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Abstract 

We consider the influence of interparticle interaction on the magnetization reversal in the oriented 
Stoner-Wohlfarth nanoparticles ensemble. To do so, we solve a kinetic equation for the relaxation 
of the overall ensemble magnetization to its equilibrium value in some effective mean field. Latter 
field consists of external magnetic field and interaction mean field proportional to the instantaneous 
value of above magnetization. We show that the interparticle interaction influences the temperature 
dependence of a coercive field. This influence manifests itself in the noticeable coercivity at T > T& 
(T;, is so-called blocking temperature). The above interaction can also lead to a formation of the 
" superferromagnetic" state with correlated directions of particle magnetic moments at T > 
This state possesses coercivity if the overall magnetization has a component directed along the 
easy axis of each particle. We have shown that the coercive field in the "superferromagnetic" 
state does not depend on measuring time. This time influences both and the temperature 
dependence of coercive field at T < T&. We corroborate our theoretical results by measurements 
on nanogranular films (CoFeB) x -(Si02)i-2; with concentration of ferromagnetic particles close, but 
below percolation threshold. 

PACS numbers: 61.18.Fs, 61.46.+W, 75.50.Tt, 75.60.Ej, 75.30.Gw, 75.75+a 
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I. INTRODUCTION 



The consideration of interparticle interaction in an ensemble of single-domain superpara- 

n n 

magnetic particles is important both from theoretical [1|, 12[ and practical points of view 
3], 0, 0, y, 13]. The above interaction alters (as compared to the case of noninteracting par- 
ticles) the magnetization curves, the coercive fields and the temperature dependence of an 
ensemble magnetic susceptibility. The intergranular interaction is always of dipole nature, 
although there are cases where an additional exchange interaction occurs also. 

To be more specific, if metallic ferromagnetic (FM) granules are embedded in conducting 
host, they can have exchange interaction via common electronic system of a composite even 
if the granules concentration is lower than the percolation threshold. In the latter case, the 
interaction is of RKKY type being small and of alternating sign. If the host is dielectric 
and the granules concentration is lower than the percolation threshold the interaction will 
be of entirely dipole nature. On the other hand, close to and/or above this threshold, the 
interaction between contacting granules may have exchange contribution of both RKKY and 
usual ion-ion types. 

A possibility for homogeneous long-range FM order to appear due to presence of dipole- 
dipole interaction has been discussed by many authors beginning from Ref. 8|. It has been 
shown in Ref. jjj, that in simple cubic lattice of magnetic dipoles the latter interaction 
generates long-range antiferromagnetic order rather then FM one. The same conclusion has 
also been drawn in Ref. Q]. The question about stability of FM long-range order in the 
cases when corresponding static mean field solution predicts the appearance of FM order 
has been considered in Ref. [lo|. This analysis shows that FM order is unstable with respect 
to 3D perturbations, but is stable with respect to 2D ones. Despite of the above discussions, 
the appearance of long-range FM order in granular systems with dielectric matrix has been 
detected experimentally in many systems with granules concentration both below and above 
percolation threshold. The properties of such FM ordered states of nanogranular systems 
have not been adequately explored. 

Regardless of the intergranular interaction nature, the question about joint influence of 
superparamagnetism and intergranular ordering on the magnetic properties of FM particles 
ensembles is still opened. It is naturally to expect that intergranular interaction generates 
a correlation of the particles magnetic moment directions. If this interaction is of FM 
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type, then below certain temperature T s f, it should generate long-range magnetic order. 
Otherwise, so-called superspinglass state can be realized. 

Most frequently, each particle has a certain crystallographic anisotropy. If its shape is 
non-spherical, the anisotropy can be of magnetostatic nature, related to demagnetization 
factors tensor for that shape. The thermal fluctuations lead to reorientation of a particle 
magnetic moment between several equivalent easy magnetization directions, dictated by the 
above anisotropy. In that case the observable magnetic properties of the ensemble are dif- 
ferent depending on the relation between the reorientation time and the period of ensemble 
observation. For noninteracting particles, the thermally activated reorientations of their 
magnetic moments are mutually independent. Two kinds of ensemble superparamagnetic 
states can be distinguished. Namely, there are equilibrium and nonequilibrium superpara- 
magnetic states. The equilibrium state occurs when the magnetic moment of an average 
particle (i.e. the typical ensemble particle with some average parameters) "covers" all per- 
missible easy magnetization directions during the time of observation. The nonequilibrium 
or "blocked" superparamagnetic state occurs in the opposite case, when the particles are 
"blocked", i.e. they cannot alter their magnetization orientations during the observation 
time. The threshold temperature between these two states is called blocking temperature, 
Tf,. The magnetic switching in the blocked state (T < T&) has a hysteretic character, while 
at T > Tft it is almost unhysteretic. 

Now we "turn on" the interparticle interaction. If it has FM character, then at T < 
T s f (T s f is determined by the interaction and can be regarded as FM phase transition 
temperature) the particles magnetic moments influence each other. Their reorientations 
cannot be independent. In such case, the particles ensemble may be considered as an 
effective ferromagnet with hysteretic magnetization reversal. The origin of hysteresis is 
similar to ordinary ferromagnets - either a pinning of the domain walls motion (if the above 
structure has domains) or simply loss of stability of homogeneously magnetized ground state 
(see, e.g. Ref. [U[ for details). 

It is naturally that can vary depending on particles size, anisotropy and observation 
time. In turn, T s f should strongly depend both on particles size and mean interparticle 
distance. Hence, the situation when either > T s f or T& < T s f can occur. The main 
problem here is possible influence of above FM ordering on temperature dependence of a 
coercive field in superparamagnetic state. Other important question is how this influence is 
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modified depending on observation time. 

The above ordered state has been identified for TJ, < T s f in Refs. Q, Q, Q. A 



nonzero, weakly temperature dependent coercive field has been observed. The authors 
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141 ] call this state "superferromagnetic" , defining it as that lying in the temperature 



range between equilibrium (unhysteretic) and blocked (hysteretic) superparamagnetic states. 
Neither temperature nor observation time dependencies of a coercive field have been analyzed 



theoretically in Refs. 



12 



13 



14|. The measurements in Refs. [12l.ll3l. Il4j were carried out on 



the films much thicker than the granules average size and with high enough relative volume 
content of the ferromagnetic component. It was supposed in Ref. [l^] that intergranular 
interaction in their films has exchange nature. At the same time, in Refs. [15I . 16], where the 



ordering in above nanoparticles ensemble has also been observed, the absence of coercivity in 
the temperature range between T s f and Tb has been reported. The coercivity appeared only 



12 



,Ll3| 



We note that in Refs. 



15 



16] 



at T < Tb in contradiction with the data of Refs. 
the materials under investigation were the multilayer films rather then above 3D systems. 
In these materials, the layers of granules with relatively small 2D filling were separated by 
the insula ting layers some thicker than the average size of a granule. It has been supposed in 



Refs. 



15 



16| that such samples emulate 2D ensembles of particles with dipole intergranular 



interaction. 

Thus, now there is a lack of complete understanding of the "superferromagnetic" state 
nature and its properties. 

In this paper, we analyse theoretically and experimentally the influence of the interac- 
tion on the coercivity of an ensemble of oriented uniaxial superpara mag netic particles. We 
consider the above interaction in a mean field approximation [l3, 111 12] without discussing 
the nature of such interaction. The magnetization switching process is considered in the 
framework of Neel discrete orientations model 20(] by the solution of a kinetic equation for 



magnetization similar to Refs. [18|, |2]J, |22j|. Contrary to those papers, here we discuss a 
possibility of superferromagnetic state creation. Also, here we study the changing of the 
temperature dependences of coercive field due to intergranular interaction with taking into 
account the influence of the observation time on the measured physical quantities. We con- 
sider the magnetization switching at linear (in time) magnetic field scanning. It corresponds 
to traditional scheme of magnetostatic measurements. 

We show that at T s f > Tb the interparticle interaction generates coercivity at T < T s f 
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and changes the ordinary (for a superparamagnet) temperature dependence of a coercive 
field at T < T&. At T < X& t he interaction yields the coercivity growth in coincidence with 
earlier Monte Carlo results 
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231 ] . We have also shown that at T < Tf, the interparticle 



interaction modifies the temperature dependence of a coercive field, obeying Neel-Brown 
law. One more result is that in "superferromagnetic" state the coercive field is related to 
the collective magnetization reversal of the particles and does not depend on measuring 
time. 



To corroborate the above theoretical results experimental 



measurements in the granular (CoFeB) x -(Si02)i-a: films 24 



y 



. we perform magnetostatic 



251 ] . In the samples under 



investigation, the FM nanoparticles were anisotropic with easy magnetization axes oriented 
along a certain direction in a film plane. We study the temperature dependences of coercive 
field at different observation times. Our experimental results are in good coincidence with 
the theoretical model presented below. 



II. THE MODEL 



Let us consider an ensemble of interacting SW-particles with their easy axes pointing 
along the direction of an external magnetic field. For this case, the energy density per 
particle in a mean field approximation is: 

U = -Kcos 2 (6) -m p (H + Am)cos(#). (1) 

Here K is a uniaxial magnetic anisotropy constant, m p is a single particle saturation mag- 
netization (it is the same for each particle), 9 is the angle between a particle magnetization 
vector and external magnetic field H direction, A is a mean field interaction parameter and 
m is an average magnetization per each ensemble particle. The latter quantity equals to the 
overall ensemble magnetization divided by the relative volume occupied by ferromagnetic 
particles in a sample. Below we will use the dimensionless magnetization M = m/m p . This 
value will be the same both for the ensemble and for each single particle. 

Let us pay attention that the potential energy profile ([1]) has the form of double-well 
potential. According to the Neel model {20I we can describe this system in the temperature 
range typical for the magnetostatic measurements (0 < T < 3Tb) as a system with two 
possible orientations of particles magnetic moments. If the magnetic field is directed along 
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easy magnetization axes of the particles, the double well potential can be substituted by its 
two lowest energy levels so that our system can be described as a two-level system. In this 
case, the transitions between levels corresponding to magnetic moment reorientations occur 
as thermally activated hops over energy barrier. 

In this model, the magnetization dynamics is of purely relaxational type. This means 
that time dependence of magnetization M at fixed temperature T and magnetic field H 



can be described (similar to Refs. 



18 



22j ) by Bloch-like equation for z-component of 



magnetization only. If the interaction term Am is absent in Eq. ([I]), the equation for 
magnetization dynamics has the form: 

^ = \ [M. - M (i) ] , (2) 

where = M(t — > oo) is the equilibrium magnetization at fixed H and T and r is the 
relaxation time. For our case of two-level system r _1 = W\ 2 + W 2 i, where Wij(j = 1,2) are 
probabilities of transition between i and j levels in the double- well potential (JTJ). According 
to approach {20! for SW particles, the final form of r reads: 

(3) 



/„ («p(-a=&) + exp(-SgS)) ' 
Here k is Boltzmann constant, /o ~ 10 8 4- 10 12 s _1 for typical magnetic particles and E\, 
E 2 , Ej, are, respectively, the energies of minima of U(9) and a barrier between them. The 
quantities E\, E 2 and E^ depend on K, A and M(t). They are the functions of time by 
virtue of M(t) dependence. The expressions fl2]) and correspond to the approach, where 
a fictitious particle (corresponding to magnetization) is localized exactly in the minimum 
of U(<?) rather then "smeared" by temperature in a wide interval of angles 9 within the 
well of the potential ([T]). For magnetic field sweeping times, taking place in magnetostatic 
measurements, this approach is well satisfied in the temperature range T < (4 4 6)71, i.e. 
in the entire temperature domain. Thus magnetization reversal occurs by thermoactivation 
overbarrier hopping. In a mean field approximation, the equilibrium magnetization for such 
two-level system is determined by usual equation: 

= tanh m p {H + \ mca )V p = tanh 2(h + X redMoD )_ 

kT T red V ; 
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Here we introduce following dimensionless parameters: the dimensionless magnetic field, 
h = H/H a (H a = 2K/m p is the anisotropy field), the dimensionless relaxation time r r = 
t/ , the temperature T red = kT/(KV p ) (V p is SW-particle volume) and the dimensionless 
parameter of interparticle interaction \ rec i = \vn? v j (2K). We also introduce the dimensionless 
energy minima Ei/(KV P ), E 2 /(KV P ), barrier maximum energy Ef,/(KV p ), dimensionless 
time t r = tfo and time of measurements t reg = t exp fr,. Here t is a real dimensional time and 
t exp is a characteristic dimensional time of measurements, i.e. "measuring time" - the time 
required for magnetic field sweeping in the range of H a . 

The introduction of the interaction term Am in the Eq.([T]) modifies the character of 
relaxation. Namely, under magnetization reversal this term becomes time dependent as it 
comprises the magnetization m(t). This means that the overall magnetization relaxes not to 
the above real equilibrium magnetization M ro = M^if, T), but to certain (so far unknown) 
self-consistent equilibrium magnetization value, dictated by the effective magnetic field H + 
Xm(t) at each time point. We denote this new hypothetical equilibrium magnetization as 
m^it), and its normalized value as M^(t r ) = m^(t r / fo) / 'm p . Here we note, that M^{t T ) 
= M^[H,T, M(t r )] so that the kinetic equation for magnetization assumes the form: 

^) = I[Mi-M W ], (2a, 
where M^(t r ) is determined by the equation 

M , Jtr) = tanh rnAH + x^/MK = tanh feMH. (4a) 

Kl J- red 



:o (Hall will be solved numerically. To model the hysteresis 



The equation (|2al) with respect 
loops, we consider, similar to Ref. 
in real experiments. In the dimensionless variables the equation (l2a|) assumes the form: 



22|, the linear field sweep h(t r ) = (t r /t reg ) — l as it realized 



dM(t r 

dt r 



exp 



(tr/Ueg ~ 2 + X re dM(t r )Y 



- red 



+ exp 



(tr/Ueg + KedM(t r )Y 



x ^tanh 



2(t r /t 



reg 



Tred 

1 + X red -M(t r ))' 



T, 



red 



- M(t r 



(5) 



The last brackets in the right-hand side of Eq. (([5])) define the difference between self- 
consistent hypothetical equilibrium magnetization ( Hal) and its current value. 
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FIG. 1: (color online) The temperature dependence of coercive field h c (y/T rec i, t reg , X re( j) for t reg = 
10 8 and different values of X red (0, 0.1, 0.2, 0.3, 0.4). Points are the results of numerical solution of 
Eq. ([5]) and the full lines correspond to low- and high temperature extrapolations of the numerical 
curves. 

III. RESULTS OF MODELING 

The solution of Eq. (jSJ) shows that the account for the intergranular interaction term 
X re dM(t r ) increases the coercivity. Also, the hysteresis loops become "more rectangular" 
("harder"). The temperature dependences of the coercive field h c , extracted from the hys- 
teresis loops calculated with the help of Eq. (jSJ), are shown on Fig. [T] at different values 
of interaction parameter X re d- It is seen that the interaction increases the coercive field. 
Besides that, at 2X re d > T b * the dependence h c (T re d) has two linear in \/T re d parts. 

The low-temperature part is similar to Neel-Brown law: 



where T fc * is a dimensionless effective blocking temperature, determined by extrapolation 
of the low-temperature part of dependence h c {\/T rei j) up to its intersection with abscissa 
axis. The value of T 6 * depends on X re d- 

The high-temperature part of h c (^T red ) lies at T red > T b *. Its origin is a consequence of 
formation (due to interaction X re d > T 6 */2) of the state with correlated directions of granules 





(6) 
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magnetic moments realized at T re d < T or d = 2X re d- Here T or d corresponds to dimensionless 
temperature of the long range magnetic ordering in the granules ensemble. In Section |T] 
(Introduction), the real (dimensional) temperature of such ordering had been denoted as 



Strictly speaking, at X re d < T£/2 and close to T red = 2\ red , the dependence h c (\/T re d) 
acquires additional slope as compared to that in Eq. ([6]). This deviation, however, cannot 
be seen in the scale of Fig. [1] so we do not plot corresponding curve on Fig. HJ 

In the temperature range above blocking temperature, the relaxation time is much less 
than measuring time so that the appearance of coercivity in this temperature range is not 
a consequence of slow system response to magnetic field sweep. The coercivity at T fc * < 
T re d < T or d is due to emergence of a self-consistent mean field of interparticle interaction 
at T re d < Tord- The direction of latter field is the same as a direction of the external 
magnetic field at the initial magnetization stage, when we lower the magnetic field from 
saturation down to zero. This mean field stabilizes the directions of particles magnetic 
moments opposite to the external field direction during field scanning process in the interval 
< h < h c - In the above temperature range T re d < T or d the overall particles magnetization 
can be well described by Eq. HJ This situation corresponds to the joint action of external 
and above self-consistent fields on each particle. In this temperature range, the quantity \h c \ 
is determined by stability limits of Eq. H] solution at opposite signs of h and X re dM. The 
values of h, where dM/dh diverges, permit to obtain the expression for the coercive field 
h c nt related to interaction term. It reads: 



Here T or d = 2\ re d is the ordering temperature. In other words, this is a temperature 
where the state with correlated directions of particle magnetizations emerges. If T ord > T 6 *, 
the above ordering occurs at T b * < T < T ord , corresponding to the parameters values shown 
on Fig. [T]and Fig. [2J The curve h c nt {T red ) determined by Eq. (j7j) has quite complex shape, 
but it is proportional to y/T re d at temperatures from 0.1T or d to 0.9T or d- In this temperature 
interval, the values of h l c nt (Eq. ([7])) can be made equal to those from Fig. [Tjat T red > T 6 *. 
Note that linear in y/T red asymptotics of h % c ni {T r ed) crosses the ordinate axis near X red . 

More detailed analysis of calculated curves presented on Fig{T] shows that they can be 
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FIG. 2: (color online) The temperature dependencies of coercive field h c (^T re d, t reg , \ re d) for 
Ked = 0.2 with different measuring times (t reg = 10 4 , 10 6 , 10 s , 10 10 , 10 12 , 10 14 respectively). Points 
and full lines are the same as on FigJH 

described with good accuracy by the sum of Neel-Brown type contribution h^ B * = 1 — 
y/T red /T^ and the contribution /i* ni determined by Eq. ([7j). Here Tj^ = T^(t reg ) is a 
blocking temperature similar to that in Neel-Brown formula but modified with respect to 
real experimental conditions (continuous magnetic field sweeping for magnetization reversal 
curves registration, see [26( for details). The value does not depend on X red . 

The dependencies h c (\/T re< i) for different measuring times t reg at fixed interaction param- 
eter \ re( i are reported on Fig. [2j One can see that the low-temperature part of the curves 
depends on measuring time owing to the corresponding dependence of parameter T^. This 
is a characteristic feature of the systems without interaction and with thermally activated 
hopping at finite measuring time. The high-temperature side of the curves does not depend 
on t reg . In this temperature region, all curves for different t reg have the same slope. 

Note that for T re( i < T ord the hysteresis curves have a "hard" (rectangular) shape. At the 
same time, for the case \ re d = and T red smaller, but close to T^, the hysteresis loops have 
"softer" or "more canted" shape. 

It is significant that the independence of a coercive field from measuring time at tem- 
peratures higher then blocking temperature is a consequence of quick (during measuring 
time or sooner) establishing of the equilibrium population in double-well potential of the 
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SW-particles in this temperature range. More detailed analysis of our model shows that 
a faint dependence h c (t reg ) is still present in this temperature range. Additionally, in this 
temperature range, the magnetization does not have a step at h = h c but rather varies 
continuously, changing its sign during the temperature dependent relaxation time r r . At 
T red > the relaxation time r r is much shorter than measuring time. Note also, that our 
model with interparticle interaction does not imply the coercivity for magnetization along a 
difficult direction. Thus, for ensemble of interacting SW-particles a temperature dependence 
of coercive field has the form: 

h c {T re d, Ked) — h^lTred, Xred) + (1 ~ \J T red /T^) (8) 

for T red < and h c (T red , X red ) = h l c nt (T red , X red ) for T red > T^, where h l c nt (T red , X red ) 
determined by Eq. (j7j). Note, that the Eq. (jSj) is exact everywhere, except for narrow region 
near where it is fulfilled approximately. 



IV. EXPERIMENTAL 



To corroborate the above theoretical approach, we measure the magnetostatic charac- 
teristics of nano-granular films (Coo.25Fe .66Bo.o9)x — (Si0 2 )i_ x . The aim was to check the 
transition of SW particles ensemble from relaxation regime of magnetization reversal to 
steady-state regime of the intergranular "ferromagnetic" ordering (arising due to the inter- 
action), when the coercive field ceases to depend on measuring time. In our measurements, 
we use the (Co .25Fe .66Bo.o9)a; — (Si0 2 )i-a; films grown in the Energy Electronics Laboratory, 
So jo University, Japan. The ferromagnetic granules were amorphous and their shape was 
close to the spherical. A strong easy-plane anisotropy related to the demagnetization fac- 
tor arose for the entire film sample. In the granular films under investigation, the uniaxial 



anisotropy in a film plane had been formed by special technological measures [24j. |25| . This 
anisotropy was supposedly related to the small deviation of the shape of granules from the 
spherical one. Thus the easiest (i.e. easy in a film plane) axes of all granules have been ori- 
ented almost parallel to each other. Therefore such granular system can be considered as an 
ensemble of easy axis oriented Stoner-Wohlfarth particles. The ensemble can be considered 
as "noninteracting" one for x substantially lower than the percolation threshold x c and as 
"strongly interacting" one for x 3> x c . According to the technologist, who had fabricated 
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FIG. 3: (color online) The temperature dependence of coercive field in the sample with x = 0.55 
measured along easy direction in a film plane. Inset reports the magnetization reversal curves 
along easy (1) and hard (2) directions in the same film plane at T = 300 K. 



our sample, the ferromagnetic component content in it is x — 0.55. It corresponds to volume 
fraction of the ferromagnetic granules about 0.26. It is strictly lower than volume fraction of 
percolation threshold. However, the film saturation magnetization, ferromagnetic resonance 
and granule magnetization 27| data have shown that real volume fraction of granules in this 
sample is essentially higher (then 0.26), up to 0.4 -r- 0.45. We assert that this real volume 
fraction is a little below percolation threshold. To prove this assertion, we had measured 
the magnetoresistance curves at T=300K. The sample has high enough specific resistance, 
p(T=300 K) = 250 mOhm/cm. The measured magnetoresistance curves contained the con- 
tribution from only tunneling magnetoresistance and did not contain the contribution from 
the anisotropic magnetoresistance. Since latter contribution in such films appears for x > x c 
only [25)], this result proves above assertion. For the above sample, the experimental data 
for magnetization in the film plane along easy direction (curve 1) at room temperature are 
presented on the inset to Fig. [3] One can see that the hysteresis loop is close to rectangu- 
lar and has almost 100% remanence. The field dependence of magnetization along a hard 
direction in a film plane (curve 2), has no hysteresis, has jogs at the intraplane anisotropy 
fields and is almost linear in magnetic field between jogs. The temperature dependence of 
coercivity (Fig. [3]) demonstrates two linear in yT parts, which coincide with the results of 
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FIG. 4: (color online) a) The temperature dependencies of coercive field H c for the sample with 
x = 0.55 versus y/T. The legend shows the field sweeping time from -500 Oe to +500 Oe and vice 
versa (t sweep ). b) The results of calculations of coercive field h c (y/T re d, t reg , X re d) f° r Ked = 0.075. 
The ratios of dimensionless measuring times t reg are the same as those on Fig. [3h. 

above theoretical modeling. 

At the same time, such dependence can be interpreted as a consequence of a bimodal size 
distribution of ensemble particles so that each particle group has its own TJ, value. To prove 
or disprove such possibility, we perform the number of magnetostatic measurements with 
different rates of magnetic field scanning. Fig. reports the temperature dependences of 
coercive field H c at temperatures from 100 to 470 K at magnetization along easy direction. 
At temperatures lower than 100 K the temperature dependence of a coercive field in this film 
demonstrates an anomaly. We will not discuss that in the present paper, having restricted 
ourselves by temperature region above 100 K only. 

The dependences H C (T) (in the form H C (VT)) on FigHk are obtained for different times 
of a magnetic field sweeping (t sweep =l min, 4 min, 16 min, 64 min) from -500 Oe to +500 
Oe and vice versa. It is seen, that increasing of measuring time (t exp is proportional to 
tsweep', actually for the anisotropy field value of this sample H a = 80 Oe, t exp = t sweep /25) 
leads, as it should be, to decreasing of H c at fixed temperature. It also yields the lowering 
of the blocking temperature TJ, determined as an intersection point of an asymptote to 
the linear low-temperature part of dependence H c (\/T) and abscissa axis. It follows from 
Fig. that the dependences H c (\/T) behave like theoretical dependences from Fig. El 
In both experimental and simulation curves we observe a noticeable coercivity practically 
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independent of measuring time above blocking temperature. The main properties of the 
presented curves can be formulated as follows. At low temperatures, T <C T& (TJ, is taken 
for the longest possible measuring time) the dependence of H c on measuring time becomes 
stronger with temperature increase. At a certain temperature, slightly lower then T&(i ea; p), 
a sensitivity of H c to measuring time variations reaches a maximum. At last, at T > TJ, 
(now Tf, is taken for the shortest possible measuring time) the H c value ceases to depend 
on measuring time. At all temperatures a hysteresis loops remain "hard", conserving a 
rectangular form. 

To compare the experimental and theoretical dependences of H c on measuring time it is 
necessary to account for the fact that measuring times, corresponding to the curves on Fig. 
121 have 10 orders of magnitude variation, while experimental data from Fig. Hh have only 64 
times difference. To illustrate the similarity between experimental and theoretical data, on 
Fig. Hb we present a number of theoretical curves with the relation of measurement times, 
identical to that in experiment. The curves are plotted for X-ed = 0.344, corresponding to 
A = 0.017, CoFeB granule magnetization m p = 1590 Gauss 23] and in-plane ahisotropy field 
H a = 2K/m p = 80 Oe. Latter value follows from the magnetization curve of our film in 
the "hard-in-plane" direction (curve 2 on inset to Fig. [3]). The value A = 0.017 is obtained 
from the equation if* n *(T — ► 0) = \m p (with respect to the value if* n *(T — > 0) = 30 Oe 
following from Fig. [3]), where if* nt is dimensional value of h 1 ™ 1 . The values of t reg (shown on 
the legend to Fig. HJd) have been chosen from the condition of best fit between model and 
experimental Tb(t reg ) values. 

Although the reported experimental results are not quite identical to the results of our 
calculation, there is obvious qualitative coincidence. Namely, in our opinion, they demon- 
strate uniquely the existence of superferromagnetic state with coercivity. This fact is also 
corroborated by Fig. [5j where the temperature dependence of relative coercivity increments 
is reported for different measuring times. We define the above increments as the difference 
of coercive fields for two substantially different measuring times divided by the coercive field 
value at larger time. This dependence is reported both for theoretical (Figj5h.) and exper- 
imental (Fig. [5b) results. It is seen, that the theoretical dependence for noninteracting 
ensemble has a sharp increase near the blocking temperature. At the same time, for the 
interacting ensemble, this dependence has a peak near Tj, with subsequent decrease. The 
experimental dependences (see Fig. |5b) also demonstrate the maximum with decrease. 
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FIG. 5: (color online) The temperature dependence of relative coercive field increments at regis- 
tration time variations: a) - theoretical results; h c i(T re( j) is taken for t reg = 10 6 and h c 2(T re d) for 
t reg = 10 8 . Full line corresponds to noninteracting case, circles - to interacting. Inset shows the 
same for wider temperature range, b) - the symbols represent experimental points, the curve is 
guide for eye; H C \(T) is taken for t sweep = 1 min and H C 2(T) for t sweep = 64 min. 

Thus it turns out that the behaviour of experimentally observed H C (T, t sweep ) dependence 
is qualitatively similar to results of our modeling for the case when intergranular interaction 
generates the state with correlated directions of particles magnetic moments and with a 
coercive field independent from measuring time. The above results illustrate the case when 
temperature T s f ("dimensional T or d") of transition to such "superferromagnetic" state ex- 
ceeds blocking temperature T&. The decreasing of intergranular interaction parameter A can 
lead to opposite situation, when T s f < TJ,. In this case the variations of dependence H c (y/T) 
also occur near T s f. However, they are too faint to be observed experimentally. 

It is possible to predict, that relaxational magnetization of the SW particles ensemble 
weakens if interaction energy exceeds the anisotropy energy. In this case the ensemble 
behaves as a uniform ferromagnetic medium with possible occurrence of "superdomains" 
consisting of many adjacent particles. 

Thus, the experimental dependence H C (T) in nanogranular magnetic film with granules 
concentration close but a bit lower than percolation threshold is in qualitative agreement 
with model predictions for the system of superparamagnetic particles with intergranular 
interaction described in a mean field approximation. Additionally, the similar results have 
been obtained in our studies of the (Coo.25Feo.66Bo.o9)z — (SiC^i-i film with nominal value 
x = 0.60. 
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Note that in spite of aforementioned qualitative resemblance of the experimental and 
model data the quantitative correspondence is not so good. The obtained t reg values on Fig. 
Sb are too small. It is the consequence of big granule size, which is needed for coincidence 
of the observed T s f and TJ, values with those expected from the model. Particularly, to 
coincide the V p value from the expression kT s f = XrripVp with A obtained from the condition 
Hl nt (T — > 0) = Xm p pa 30 Oe, we need the mean diameter of granule 14 nm. To coincide 
the observed values of T& and H a under usual assumption \og(t exp fo) = 20-r25, we need this 
diameter to be about 24 nm. These mean diameter values are essentially more than those 
expected in film fabrication process. One of possible explanations of such discrepancy is the 
relative simplicity of used theoretical models. 

The last (but not the least) question is about the nature of intergranular interaction in 
the studied samples. Is it of dipole-dipole or exchange nature? We do not have a convincing 
answer to this question. One of possible suppositions is that exchange part of the interaction 
coexists with a dipole-dipole one due to closeness of granules content in our film to the 
percolation threshold. 

V. CONCLUSIONS 

To conclude, here we present a mean-field consideration of the magnetization of ensemble 
of interacting Stoner-Wohlfarth particles. We do that on the base of the kinetic equation 
solution. The equation has been written for the relaxation of overall ensemble magnetization 
to its self-consistent equilibrium state in the effective field consisting of external and the 
interaction fields. The latter field, in turn, is proportional to instantaneous value of overall 
magnetization. Numerical solution of the above kinetic equation shows that the presence of 
mean-field interparticle interaction leads to the following effects: 

- At certain temperature, T s f, proportional to interaction parameter A, the system of 
FM granules undergoes the intergranular magnetic ordering - "superferromagnetism", yield- 
ing the additional coercivity at T < T s f. For T s f > T b , the essential coercivity arises at 
temperatures above blocking temperature; 

- At T < T s f a coercive field increases and magnetization reversal becomes "harder", i.e. 
the hysteresis loops become almost rectangular with increased remanence; 

- Temperature dependence of coercive field in the low temperature region resembles very 
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much Neel-Broun law; for T s f > Tb at increasing temperature the dependence H c (yT) (or 
h c (y/T Tec i) in the dimensionless units) has an inflexion point at blocking temperature and 
then continues to go linearly in \/T up to T s f with much smaller slope; 

- In the system of weakly interacting particles the values of TJ, and H C (T) at T < T& 
depend on the measuring time as it is usual for SW particles. At the same time, in the 
case of T s f > Tb the dependence H c (\/T) in the range Tb < T < T s f ceases to depend on 
measuring time. 

- The temperature dependence of a coercive field (Eq. ([8])) is described with good accu- 
racy by two additive contributions. The first one (below blocking temperature) is strongly 
dependent on measuring time and reflects coercivity related to a metastability of the system 
at finite measuring times. The second one, which is almost independent from measuring 
time, has a temperature dependence described by Eq. (J7j). This contribution reflects a 
change in a mean field of intergranular interaction in the process of magnetization reversal. 

All above manifestations of interparticle interaction in SW particles ensemble are ob- 
served experimentally in magnetostatic (with 64 times difference in measuring times) 
measurements of a magnetic field and temperature dependencies of magnetization of 
(Coo.25Fe .66Bo.o9)o.55 — (Si0 2 )o.45 nanogranular films with FM granules content close, but 
below a percolation threshold. Thus it is firmly established that within the described ap- 
proach the results of our numerical modeling, are in good coincidence with the experimental 
data. 

Acknowledgments 

This work was partly supported by the grant of NAS of Ukraine Target Program " Nanos- 
tructural systems, nanomaterials and nanotechnologies" . 



[1] Z. Mao, D. Chen, and Z. He, J. Magn. Magn. Mater. 320, 642 (2008). 
[2] A. D. Liu and H. N. Bertram, J. Appl. Phys. 89, 2861 (2001). 

[3] C. Binns, M. J. Maher, D. Kechrakos, and K. N. Trohidou, Phys. Rev. B 66, 184413 (2002). 
[4] P. Allia, M. Coisson, P. Tiberto, F. Vinai, M. Knobel, M. A. Novak, and W. C. Nunes, Phys. 
Rev. B 64, 144420 (2001). 

18 



[5] J. Escrig, S. Allende, D. Altbir, and M. Bahiana, Appl. Phys. Lett. 93, 023101 (2008). 
[6] M. Hillenkamp, G. Domenicantonio, and C. Felix, Phys. Rev. B 77, 014422 (2008). 
[7] D. Yao, S. Ge, X. Zhou, and H. Zuo, J. Appl. Phys. 104, 013902 (2008). 
[8] J. M. Luttinger and L. Tisza, Phys. Rev. 70, 954 (1946). 
[9] E. Z. Meilikhov and R. M. Farzetdinova, JETP 94, 751 (2002). 
[10] Y. G. Pogorelov, G. N. Kakazei, M. D. Costa, and J. B. Sousa, J. Appl. Phys. 103, 07B723 
(2008). 

[11] S.V.Vonsovsky, Magnetizm (John Wiley, New York, 1974). 

[12] R. S. Iskhakov, G. I. Frolov, V. S. Zhigalov, and D. J. Procof'ev, Lett. J. Tech. Phys. 30, 51 
(2004). 

[13] B. J. Jonsson, T. Turkki, V. Strom, M. S. El-Shall, and K. V. Rao, J. Appl. Phys. 79, 5063 
(1996). 

[14] J. P. Perez, V. Dupuis, J. Tuaillon, A. Perez, V. Paillard, P. Melinon, M. Treilleux, B. Barbara, 
L. Thomas, and B. Boushet-Fabrer, J. Magn. Magn. Mater. 145, 74 (1995). 

[15] W. Kleemann, O. Petracic, C. Binek, G. N. Kakazei, Y. G. Pogorelov, J. B. Sousa, S. Cardoso, 
and P. P. Freitas, Phys. Rev. B 63, 134423 (2001). 

[16] X. Chen, S. Sahoo, W. Kleemann, S. Cardoso, and P. P. Freitas, Phys. Rev. B 70, 172411 
(2004). 

[17] D. L. Atherton, IEEE Trans. Magn. 26, 3059 (1990). 

[18] M. A. Chuev and J. Hesse, J. Phys.: Cond. Matt. 19, 506201 (2007). 

[19] J. J. Zhong, J. G. Zhu, Y. G. Guo, and Z. W. Lin, IEEE Trans. Magn. 41, 1496 (2005). 

[20] L. Neel, Ann. Geophys. 5, 99 (1949). 

[21] M. A. Chuev, JETP Lett. 85, 744 (2007). 

[22] A. A. Timofeev, V. M. Kalita, and S. M. Ryabchenko, Low Temp. Phys. 34, 446 (2008). 
[23] L. Wang, J. Ding, H. Z. Kong, Y. Li, and Y. P. Feng, Phys. Rev. B 64, 214410 (2001). 
[24] M. Munakata, M. Yagi, and Y. Shimada, IEEE Trans. Magn. 35, 3430 (1999). 
[25] P. Johnsson, S. I. Aoqui, A. M. Grishin, and M. Munakata, J. Appl. Phys. 93, 8101 (2003). 
[26] A. A. Timopheev and S. M. Ryabchenko, Ukr. J. Phys. 53, 261 (2008). 

[27] M. Munakata, M. Namikawa, M. Motoyama, M. Yagi, Y. Shimada, M. Yamaguchi, and 
K. Arai, J. Magn. Soc. Japan 26, 388 (2002). 



19 



